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Summary 


The onset of electrohydrodynamic motion associated with the imposition of an electric field 
across a thin layer of liquid has been investigated for the case in which the electrical 
conductivity varies linearly over the depth of the layer. The variation of the conductivity is 
due to concentration gradients in the charge-carrying solutes and its spatio-temporal evolution 
is represented by a convective-diffusion equation. When the viscous relaxation time is long 
compared to the time for charge relaxation, the analysis reveals that the neutral stability 
curves for the layer can be characterized by three dimensionless parameters: 
Ra e = deE*A(j/nK cff <r 0 , an electrical Rayleigh number; A<t/<t 0 , the relative conductivity 

increment, and oc, the transverse wave number of the disturbance. Here d is the thickness, 
£ is the dielectric constant, and fl is the viscosity of the layer, E 0 is the applied field 
strength at the lower conductivity boundary, and K el( is an effective diffusivity associated 
with the Brownian motion of the charge-carrying solutes. With viscous-stress-free 
boundaries, at which the electrical conductivity and the normal component of the electric field 
are prescribed, the critical Ra e is 1.504 x 10 4 at a critical transverse wave number of 1.98 
when A<t/<t 0 is 10. As Ao/Oq increases, the critical Ra e increases and shifts to shorter 
wavelength disturbances; the critical imposed field strength, however, passes through a 
minimum because the lower-conductivity boundary exerts a considerable stabilizing influence 
in the presence of steep conductivity gradients. Similar trends were obtained for liquid layers 
with rigid boundaries. 
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Chapter 1 : Background and Motivation 


Background 

Developments in the area of biotechnology have created the need to isolate 
and purify large quantities of biological materials. A common (analytical) 
technique used to purify small quantities of biomaterials is electrophoresis. 


Electrophoresis is the single most useful technique for the analysis and 
separation of proteins and nucleic acids. Its usefulness is based on its high 
resolution and applicability to a broad range of entities, from simple ions to 
living cells. It is also an inherently gentle technique, without deleterious 
effects on such delicate biological structures as organelles, living cells or intact 
chromosomes f 1 ]. 


Electrophoresis is based on the principle that most biological materials in 
an aqueous solution are electrically charged (due to ionization) and, 
therefore, move under the influence of an external electric field. 
Biological particulates and solutes migrate at a rate that is dependent on 
such things as their electrical charge, their size and shape, and the 
properties of the solution in which they are immersed. Biomaterials for 
which electrophoresis may be employed include macro-ions, complex 
polyelectrolytes, and colloidal particles and living cells. Analytic 
electrophoretic separations are carried out in four classical modes: moving 
boundary electrophoresis (MBE), zone electrophoresis (ZE), 
isotachophoresis (ITP), and isoelectric focusing (IEF) [2], 

With MBE, ZE and ITP, isolation and separation of species is based 
on differences in the electrophoretic mobility of the components — that is, 
on the rate of migration of the species in solution under the influence of an 
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external electric field. In IEF, on the other hand, the separation is based 
on differences in the isoelectric points of components [3], IEF was 
developed to fractionate amphoteric species, molecules that act either as 
acids or as bases depending on the local pH of the solution. Amphoteric 
species typically contain both carboxyl and amine groups, so IEF is used on 
mixtures of peptides and/or proteins. 

IEF depends on the formation of a stable pH gradient along which 
the components migrate until reaching their respective isoelectric points 
(pi). The pi of a charged species is the pH at which the net surface charge 
is nil [2]. An externally-applied electric field acts on the charged 
component and displaces it toward the region of isoelectric pH. Any 
movement away from this point in the pH gradient results in the ionization 
of the species, which causes the species to migrate back to its pi. Thus 
steady-state is achieved. 

Though IEF is a technique tailored for the separation and 
characterization of amphoteric materials, when carried out in free solution, 
its effectiveness is diminished by fluid motion that disrupts the pH gradient. 
Experiments were conducted at The University of Arizona Center for 
Separation Science (CSS) to examine the causes of convection in IEF. The 
goal of The University of Arizona work was to develop ways to suppress 
the deleterious effects of convection so that free-solution IEF could be 
scaled-up to the preparative level, permitting large quantities of peptides 
and proteins to be processed. The CSS experiments showed that, owing to 
concentration gradients, buoyancy-driven convection occurred when the 
solutes to be separated were localized in considerable quantity — a problem 
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termed zone over-loading. Sedimentation occurred if the biological 
particulates were too large. Furthermore, imposition of an electric field 
caused joule heating in the electrolyte and a radial temperature gradient 
developed, causing natural convection. Joule heating limited the cross- 
sectional area of the device in which a quiescent IEF separation could be 
carried out to that of a capillary. Obviously such device-size limitations 
meant that free-solution IEF could not be done on a preparative scale 
unless measures were taken to either control or circumvent the problems 
associated with natural convection. 

Subsequent investigations of IEF by the CSS (under the auspices of 
NASA) made use of microgravity conditions to suppress buoyancy-driven 
flow. The microgravity experiments uncovered still another form of 
convection that disrupted the focusing process. It was postulated that the 
convection observed in the microgravity experiments was due to 
electroosmotic and electrohydrodynamic effects. These electrically-driven 
flows, though present in the terrestrial experiments, were masked by 
gravitational effects. 

Electroosmosis arises as a result of the electrical charge present on 
the surfaces of the IEF chamber. Charge on the glass vessel walls interacts 
with the ions in solution, and this causes a diffuse layer of charge to 
accumulate immediately adjacent to the chamber walls. The externally- 
applied electric field acts on the accumulated charge, setting the electrolyte 
in motion. Electrohydrodynamic effects are caused by local deviations 
from electroneutrality within the volume of the solution. Free volume 
charge accumulates due to conductivity gradients that develop during IEF. 
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The presence of conductivity gradients is a characteristic feature of IEF 
and can be understood qualitatively by recognizing that zones in which 
different species concentrate will conduct current at different rates. As 
shown in Fig. 1 , quantitative descriptions of the conductivity profile in an 
IEF chamber are available through the use of a simulation program 
developed at the CSS. 

One clear implication of the CSS experiments is that even if 
buoyancy-driven convection can be eliminated — either by use of space- 
based platforms, or by innovative instrument design — free-solution IEF 
cannot be done on the preparative scale unless electrically-driven 
convection is also controlled. A more subtle implication of the CSS 
experiments is that producing quiescent solutions for preparative-scale IEF 
may be inherently problematic. Indeed preparative-scale IEF instruments 
developed at the CSS make use of forced convection to over-ride natural 
and electrically-driven convection; that is, preparative-scale IEF was shown 
to be quite feasible if the focusing process was carried out in solutions that 
were made to flow by design rather than by deleterious effects. 

As noted, steep electrical conductivity gradients develop along an 
IEF chamber as the focusing process proceeds. Detailed analysis of the 
flight data suggests that such conductivity gradients are the primary cause 
of the undesirable convection observed in the space experiments. A 
question that ought to be addressed, then, is the following: Suppose that a 
perfectly quiescent solution can be produced for the initial stages of an IEF 
experiment. As the conductivity gradients in the solution continue to 
sharpen, will the solution remain quiescent? It is the purpose here to 
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a) b) c) 





COLUMN LENGTH (cm) 


Fig. 1 Computer simulated concentration (a), pH (b) and conductivity 
profiles for a mixture of Arginine (Arg), Histidine (His) and 
Glutamic acid. Time points are, from bottom to top, 0, 5, 10, 
15, 20, 50, 100, 150 and 400 min. Offset for concentration, 
pH and conductivity are 50 mM, 3.4 pH units and 0,07 S/m 
respectively [3]. 
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Table 1 

Computer Simulation 


General Values 


Column length 
Segmentation 
Current density 

1 cm 

100 segments/cm 
0.001 A/cm ^ 



pK a and Mobility Values 

Component 

pK a l 

pKa2 

Mobility x 10 4 
[cm 2 /Vs] 

Glutamic acid 

2.16 

4.29 

2.97 

Histidine 

6.02 

9.17 

2.85 

Arginine 

9.04 

12.48 

2.26 

H 30 + 



36.27 

OH' 



19.87 
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examine this very question by investigating the onset of 
electrohydrodynamic motion associated with sharp conductivity gradients. 

For the balance of this chapter, the results of the CSS space-flight 
experiments will be discussed, beginning with a summary of how an IEF 
experiment proceeds. Next, the flight results will be recounted and, 
thereafter, results from computer simulations of the flight experiments will 
be presented. Having thus established a feel for the physico-chemical 
process involved in IEF, we then proceed to Chapter 2, in which an 
analysis of the onset of electrohydrodynamic motion is performed. Note 
that data and figures pertaining to the IEF space flight experiments have 
been reproduced from articles and technical reports written at the CSS [4], 

Analysis of Space Flight Experiments 

Summary Description of IEF 

The general mechanism of IEF may be described with the help of Fig. 2, 
where the separation mechanism for a three component mixture is shown. 
The most basic species of the mixture is component A, the most acidic is 
component C, and component B has an intermediate pi. 

Two separate phases may be identified during the process: a 
separation phase and a stabilizing phase. In the first phase, components in 
solution distribute in zones along the capillary separation length according 

to their pis. In the second phase, the system stabilizes and reaches steady- 
state. 

During the separation phase, two distinct transient stages are seen. 
Fig. 2a depicts the initial distribution of components in the separation 
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capillary, before any external electric field is applied; note the components 
are distributed evenly. When the potential is applied (Fig. 2b), two moving 
boundaries are formed near each electrode. An a boundary appears 
between pure-component and mixed-component zones. Faster moving (5 
boundaries appear between mixed component zones. The arrows beneath 
the boundaries indicate the direction of motion of the species; the length of 
the arrow specifies the relative speed of boundary movement. Thus, A and 
C are the fastest moving components, and they move toward opposite 
electrodes. In Fig. 2c, the faster ft boundaries meet, and the original three 
component mixture has vanished. A pure B component zone then appears 
(Fig. 2d) and again there are four moving boundaries. The outer (a) 
boundaries migrate inward, reflecting motion of component B toward its 
pi. The inner (y) boundaries move toward the electrodes, indicating that 
components A and C tend to focus at the electrodes. As a and y 
boundaries meet, a three zone solution forms, with each zone comprised of 
a pure component (Fig. 2e). In this stage, 8 boundaries form and migrate 
slowly toward the electrodes. This marks the end of the separation phase 
and the beginning of the stabilizing stage. During the interval between 
stage (e) and (f), the steady-state composition distribution is slowly 
approached. The stabilizing phase is characteristically a much slower 
process that the separation phase; typically the same time is required to pass 
from stage (a) to (e) as is then required to go from stage (e) to (f) [3]. 



Flight Experiments 

IEF experiment were performed in microgravity conditions during two 
NASA missions: STS-1 1 and STS-26. 

The first space experiment was aimed at exploring the influence of 
electroosmosis on IEF. Electroosmotic effects were studied by applying an 
anti-electroosmosis coating on the inside of the chamber walls, and by 
segmenting the separation chambers so as to interrupt flow along the walls. 
IEF chambers were built from segments of glass conduit, which were 
attached end to end in three ways: (1), as a simple series of cylinders; (2), 
as a series of tubes connected by monofilament screens; and (3), as a series 
of tubes with mylar constrictors set in between (i.e., baffles of rigid mylar 
with openings half the internal diameter of the glass segments). The five 
chamber types are shown in Fig. 3. The current passing through the 
chambers was recorded to monitor conductivity changes in the solution. 
Also, photographic records of the protein focusing were obtained since 
"red" (Hemoglobin) and "blue" (Albumin) stained species were chosen as 
the compounds for separation. 

Plots of normalized current versus time follow in Figs. 4-9. The 
experimental results indicate that the columns partitioned by monofilament 
screens perform best. With monofilament screens the current decreases 
steadily as the charged species reach their pis, and the final current is about 
one tenth of the initial value. The columns with mylar baffles also 
exhibited focusing that occurred steadily and completely. However, 
columns that did not contain partitions displayed a quite different behavior: 
a sudden increase in current followed an initial decline, and the final 
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Fig. 3 Cell cross-sections for the STS-1 1 experiments. 
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Fig. 4 Plot of current measurements versus time 
for the cells on the STS-1 1 mission. 
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Fig. 5 Plot of normalized current measurements versus 
time for the cells on the STS-1 1 mission. 


% of Initial Current 



Amonoime hocusing Solutions 

Cell ' 2 : coated 

Cell - I - uncoated 

Cell * 5 : segmented 

Cell s 5 compartmented with mylar 

ceil - 4 ■ compartmented with screens 


Fig. 6 Plot of normalized current measurements versus time 
for the cells containing solutions of ampholine. 
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Fig. 7 Plot of normalized current measurements versus time 

for the cells containing solutions of 3-component buffer. 
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Fig. 8 Plot of normalized current measurements versus time 
for non-compartmented cells on the STS- 1 1 mission. 
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Fig. 9 Plot of normalized current measurements versus time 
for segmented cells on the STS-1 1 mission. 



current measured close to the initial value. The same occurred when anti- 
electroosmosis coatings were present. Since the current increases were not 
synchronized for these columns (Figs. 8 and 9), an external disturbance 
such as a sudden acceleration of the Shuttle does not explain the data. 
Rather the onset of electrically-driven convection within the separation 
chamber, causing mixing of the components, seems to be the most plausible 
explanation for the current increases. 

Although the chambers with baffles or screens performed well, the 
fluid velocity observed in the unsegmented chambers was ca. 1 cm/min., 
much too strong to be attributed solely to electroosmosis. Thus, the 
current increases must have involved another kind of electrically-driven 
fluid motion, which Rhodes et. al. [5] subsequently deduced must be 
electrohydrodynamic in nature. 

Another noteworthy observation was that once convection ensued, it 
continued throughout the experiment; a more intuitive, cyclic behavior of 
focusing-mixing-refocusing did not occur. 

Once all elements from the STS-1 1 experiments and their post-flight 
analysis were put together, a new set of experiments was devised and flown 
on STS-26. The STS-26 experiments were aimed at verifying the STS- 11 
results, and at defining the relative magnitude of electroosmotic (EO) and 
electrohydrodynamic (EHD) effects. EO is a wall effect (proportional to 
the surface area of the IEF wall) that is proportional to the magnitude of 
the applied field. On the other hand, EHD is a bulk effect (proportional to 
the volume of the column), and it is proportional to the square of the 
magnitude of the applied field. Thus, varying the surface to volume ratio 



of the IEF chamber at constant field strength varies the relative strength of 
the two flows. The chamber configurations tested on STS-26 are shown in 

Fig. 10. Note that cylindrical chambers were included to repeat the STS- 
11 experiments. 

The experimental data from STS-26 are summarized in Figs. 11-16. 
The extent of the focusing, the time of first observable convection, and the 
strength of flow were comparable to the STS- 11 observations for the 
cylindrical chambers. Screens and a larger surface to volume ratio 
improved the focusing performance of the columns, implying that EHD is 
the principal cause of the STS- 11 current increases. Moreover, the 

performance of the rectangular cells showed a surface to volume ratio 
dependence. 


Simulation Results 

The simulation package developed at the CSS was used to estimate the 
magnitude of the conductivity gradients that developed during the STS- 11 

and STS-26 flight experiments. Refer to Table 2a for data used in the 
simulation. 

The evolution of conductivity profile obtained is depicted in Fig. 17, 
and the final concentration profile in Fig. 18. The results obtained show 
evidence of numerical instabilities in the computer code and the simulation 
failed before reaching steady-state. Problems arise in the simulation 
because of the steep pH gradients that develop. Also, the complexity of the 
separation of a five component mixture adds to the numerical difficulties. 
In order to obtain a better estimate of the conductivity profile, it was 
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Cell cross-section for the STS-26 experiments. 
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Plot of current measurements versus time for 
cylindrical cells on the STS-26 mission. 
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Fig. 12 


Plot of normalized current measurements versus time 
for cylindrical cells on the STS-26 mission. 
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Fig. 1 3 Plot of normalized current measurements versus time 
for cylindrical cells on the STS-26 mission. 
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Plot of current measurements versus time for 
rectangular cells on the STS-26 mission. 
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Fig. 16 


Plot of normalized current measurements versus time 
for cylindrical cells on the STS-26 mission. 
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Table 2a 

Computer Simulation Data for Flight Experiments on STS-11 


General Values 


Column length 
Segmentation 
Current density 

4.5 cm 
100 segments 
336.7 A/m 2 * 




pK a and. Mobility Values 




Mobility x 10 8 

Cone, x 10 3 

Component 

pK a l 

pKa2 

[m 2 /Vs 1 

[Ml 

Arginine 

9.04 

12.48 

2.71 

8.0 

Lysisl-Aspartate 

4.1 

7.7 

2 

10.0 

p-ABA 

2.41 

4.85 

3.28 

6.0 

Hemoglobin 




0.0382 

Albumin 




0.0110 

* This implies a constant voltage of 75.0 V. 


Conductivity, a (S/m) 


u.w 

0.03 

0.02 

0.01 

0 



Column length, z (m) 



Ampholytes Concentration (M) 
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Fig. 18 Simulated concentration profile at 6 min into the simulation for 

the system in Table 2a. The simulation failed after the sixth minute. 
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necessary to re-run the simulation under conditions of reduced voltage. 
Furthermore, because small amounts of protein do not significantly alter 
the conductivity profile, the revised simulation was done for a three 
ampholyte solution {sans proteins). Input data for the revised simulation 
are given in Table 2b and the results are shown in Figs. 19-23. The 
evolution of the conductivity profile is depicted in Figs. 19 and 20. Here 
the rate of focusing of the components is slower than in the flight 
experiments because of the smaller current used. Notice that it is possible 
to distinguish the migrating zones at the different stages of the process. 
More importantly, recognize that at some time between 3 and 4 minutes 
into the simulation, a large conductivity gradient develops as the pure 
Lysil- Aspartate zone appears. These localized conductivity changes are 
accompanied by free charge accumulations, and it is at this stage that the 
conductivity gradients are postulated to favor the onset of EHD. The 
magnitude of the normalized conductivity gradient, Vcr/cr, is greater than 
100 nr 1 , as shown in Fig. 21. An estimate of the minimum conductivity 
gradient required to yield flow strengths of the order of those observed in 
the flight experiments is Vcr/cr « 3 w _1 (cf. Appendix 1). The simulations 
clearly show that, if the IEF solution were to remain quiescent, values for 
Vcr/cr far in excess of 3 would develop (cf. Fig. 21). Thus, the IEF process 
proceeds in such a way that the requisite conditions are manifested for 
strong EHD floor. In Chapter 2, therefore, we take up the question of the 
onset of EHD convection in liquids that contain appreciable conductivity 
gradients. Note for completeness the steady-state pH, conductivity and 
concentration profiles are shown in Figs. 22 and 23. 


Table 2b 

Cojnguter_S^ for Flight Experiments on STS-11 


General Values 


Column length 
Segmentation 
Current density 

1.0 cm 

200 segments 
20 A/m 2 * 




pK a and Mobility Values 

Component 

pK a l 

pK a 2 

Mobility x 10 8 
[m 2 /Vs] 

Cone, x 1(P 

m 

Arginine 

9.04 

12.48 

2.71 

8.0 

Lysisl-Aspartate 

4.1 

7.7 

2 

10.0 

p-ABA 

2.41 

4.85 

3.28 

6.0 


* This implies a constant voltage of 11.53 V. 
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CHAPTER 2t : ELECTROHYDRODYNAMIC INSTABILITY IN 

a Thin fluid layer with an 
Electrical Conductivity Gradient 


Introduction 

Electrohydrodynamic motion is known to occur in liquids where spatial gradients in the 
electrical properties are present. Classical examples of these flows are found in systems 
with step-changes in the electrical properties, such as those at the boundary between 
immiscible phases. Taylor and McEwan 1 analyzed the stability of a horizontal interface 
between two superposed liquids, and Michael and O'Neill 2 investigated the stability of a 
non-conducting layer sandwiched horizontally between two conducting fluids, in a vertical 
electric field. These studies show that, without the stabilizing effect of interfacial tension, 
the quiescent fluid arrangement is disrupted. Indeed, in a series of experiments on miscible 
laminar jets, Rhodes, et al } demonstrated that, unimpeded by the restoring influence of 
interfacial tension, electrohydrodynamic stresses flatten cylindrical jets into ribbons. Still 

more examples of flows driven by interfacial electric stresses are reviewed by Melcher and 
Taylor 4 . 

Electrohydrodynamic motion has also been studied in systems where the spatial 
variations in the electrical properties are far less dramatic. Turnbull and Melcher 5 examined 
the stability of a perfectly insulating fluid layer in the presence of thermally-induced 
gradients in density, dielectric constant and viscosity. Variations in the dielectric constant 
give rise to dielectrophoretic body forces on the fluid and free charge accumulates in the 
bulk. It was shown that if a dc electric field is applied to the (quiescent) fluid layer, 

t Extensive excerpts of this chapter have been submitted in June 1994 to Physics of Fluids for publication. 


electrohydrodynamic instability sets in, provided that the thermal diffusion time is long 
compared with the charge relaxation time. Melcher and Firebaugh 6 studied fluid motion in 
semi -conducting liquids with thermally-induced conductivity gradients. They showed that, 
subject to a potential wave that travels orthogonal to the temperature gradient, the liquid 
flows anti-parallel to the conductivity gradient. 

Here we consider the onset of electrohydrodynamic motion in a layer of liquid 
across which there is a continuous variation in the electrical conductivity; in contrast to 
Melcher and Firebaugh 6 , we are concerned with conductivity variations associated with 
gradients in the concentration of charge-carrying solutes. A linear stability analysis is 
performed on the quiescent base state that exists when the conductivity varies linearly with 
position and there is an electric field applied across the layer of liquid. The conditions for 
marginal stability are determined, and the analysis shows that the relevant dimensionless 
groups are: Ra c , an electric Rayleigh number; Ao/tx 0 , the relative conductivity increment; 

r, the ratio of viscous to charge relaxation times; and oc, the transverse wave number of 
the disturbance. Numerical calculations show that, for a given relative conductivity 
increment, the fluid layer becomes unstable when the electric Rayleigh number exceeds a 
certain value; the critical values for Ra e and a depend solely on Act/ct 0 in the limit 

T °°’ which are the conditions for ohmic, aqueous systems. 


The balance laws include the effects of diffusive processes associated with the 
conductivity profile and, so, the analysis differs from those of Hoburg and Melcher 7 and 
HoburgS for the stability of fluid layers with colinear field and conductivity gradients. In 
their analyses, the conductivity gradients developed (quiescently) due to the diffusion of 
charge-carrying solutes, but these diffusive processes were then omitted from the equations 
governing the stability of the fluid layer. The rationale for such an approximation is that the 
time scale for diffusion of ionic solutes is typically quite long compared to the vis 
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response time, i.e. K ei( /v« 1, where K e{{ is an effective diffusivity for the conductivity 
[cf. Eq. (12)] and v is the kinematic viscosity of the liquid. 

In the present context, however, the electrical conductivity gradients are large. 
Thus, diffusive processes are important in spite of the fact that K ef( /v«l. More 

specifically, diffusion of the ionic solutes tends to stabilizes the layer against the onset 
electrohydrodynamic motion. To see this, consider the fluid layer depicted in Fig. 1, 
where the electric field is directed parallel to the conductivity gradient. If a parcel of fluid 
within the layer is displaced upward, it moves to a region in which, locally, the 
conductivity is higher and the electric field is lower. It can be shown 9 that a low 
conductivity parcel surrounded by higher conductivity fluid will tend to migrate toward 
regions of lower field strengths. Such migration is driven by the gradient in the electric 
field and is referred to here as dielectrophoretic, though the effects are due to variations in 
the conductivity and not the dielectric constant of the medium. These dielectrophoretic 
effects favor continued upward movement of the fluid parcel. On the other hand, diffusive 
transport of charge-canying solutes tends to eliminate the conductivity differences between 
the upwardly directed fluid parcel and its surroundings. If this diffusive transport is 
comparatively rapid, the upward motion of the parcel is arrested because the 
dielectrophoretic forces vanish with the conductivity difference. Note the fluid viscosity 
also tends to resist continued motion of the parcel. Depending on the relative strengths of 
these competing effects, which are embodied in Ra e , the displaced fluid parcel may or may 
not continue upward— and this distinguishes unstable from stable base states. If K e(( /v is 
set to zero, as was done by Hoburg and Melcher? and Hoburg,** the stabilizing influence of 
the diffusive processes is not captured and the fluid layer is made linearly unstable by the 
imposition of the external field. Such a result is analogous to what one would expect from 
the classical Rayleigh-Benard problem if, in writing the thermal energy balance for the 
perturbation variables, one were to consider the fluid to have zero thermal conductivity. 
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Interest in the present problem stems from electrically-driven convection observed 
in space-flight experiments on isoelectric focusing 10 , a separation process used to isolate 
and purify biological materials 11 . Experimental measurements and numerical 
simulations 10 12 show that appreciable conductivity gradients develop during isoelectric 
focusing (IEF), and it has been postulated that these gradients precipitate the observed 
flows ’ . The conditions leading to the onset of motion in IEF are of concern because 
convection disrupts the separation process; the linear stability analysis is formulated here so 

as to describe conditions that give rise to such deleterious convection in initially quiescent 
fluid layers. 


The presentation is organized as follows. First, the balance laws that govern the 
electrohydrodynamics of the fluid layer are presented, and then, after a brief scale analysis, 
they are placed in dimensionless form. Next, a linear stability analysis is performed; the 
quiescent base state is perturbed slightly, and the perturbations are examined. Finally the 
results of the analysis are discussed and concluding remarks are made. 


Balance Laws 

The electrical force acting on ionic solutes is transferred to the solvent through collisions at 
the molecular level. The free charges individually transfer their momentum to the solvent, 
so the Navier-Stokes equations, which govern the fluid motion, must be modified to read 

P (^ + V ' Vv ] = “ V/? + //V2v + f e ’ (I) 


where f e is the electrical force on the solvent per unit volume. Provided that the electric 
displacement D and the electric field E are co-linear, f e can be expressed in terms of the 
Maxwell stress tensor 13,14 , viz. 
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( 2 ) 

( 3 ) 

f e =PfE-i(EE)V £ . (4) 

The first term on the right hand side of Eq. (4) is the (coulombic) electrical body force 
caused by the action of the electric field on p { , the free charge density in solution. The 

second term accounts for electric forces due to gradients of the dielectric constant £ . When 
Eq. (4) is combined with Eq. (1), the momentum balance becomes 

P-^ = - v P+P,E-i(E.E)V£ +#1 v V (5) 

with V • v = 0 . (6) 

Analysis of Eq. (5) is often difficult because, according to Maxwell's equations, the 
electric field is coupled to the free charge density p f . Moreover, the free charge density is 
influenced by convection, so one may not generally specify p f and E independently of v . 
As long as the electrical current density J f is modest, E is quasi-static, so that 

V x E = 0 (7) 


f =V T, 


where for an incompressible linear dielectric, 


T = ED-^-(D-D) I . 


Since D = £E and p f = V • D, 


and 
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V-(flE)=p f . 


( 8 ) 


This latter relation is Gauss' law. Conservation of charge requires 

^ L + V-J f = 0. 

dt 1 

For ohmic conductors the current density is given by 13,14 


Jf-crE + vp f , (io) 

where <T is the electrical conductivity and the second tern, on the right hand side represents 
convection of free charge. Combining Eqs. (9) and (10) yields 

Dp f _ 

= Va • E + <r V ■ E , nn 


since V • v = 0 . 

Because the solution is an ionic conductor, the conductivity depends on the local 
ion concentration. As a result, the conductivity may vary with position and time. 

Melcher has shown that the evolution of the conductivity profile is satisfactorily 
described by 


Dg 

~Dt 


= K e{{ V 2 cr, 


( 12 ) 


wnere rt e(f represents an effective diffustvity that arises due to Brownian motion of the 

Though Eq. (12) appears to be a balance on a physical property it should be 
recognized that, in ionic conductors the conductivity varies locally with the concentration of 
charge-carrying solutes. Equation (12) holds'* provided that the local charge 
accumulations relax raptdly compared to the time necessary for momentum transfer 
(viscous relaxation time) and the time required for ion electro-migration, viz. 


£ 

« 




d d 2 

~ « — and 
(oE n v 


d 2 

Q)k B T ’ 


where tu is a characteristic mobility of the charge-carrying solutes and V is the 
Boltzmann temperature. 

Scale Analysis 

As an IEF experiment proceeds, concentration, pH and conductivity gradients develop 
along the axis of the imposed electric field. Simulations show 1 '.'2.15 [hat _ jf ^ process 

proceeds without convection, the local conductivity can vary by more than an order of 
magnitude over a length of 1mm . The dtameter of a cylindrical IEF separation chamber is 
roughly 1cm and the chamber is several cm long, so the axial distance over which the 
conductivity gradients occur is thin compared to radial and axial dimensions of the 
chamber, and the region of variation is distant from either electrode 1 1 - ,: 2 - 15 

For the situation depicted in Fig. 1 , the conductivity profile <7 V) varies linearly 
across the domain 0 < z <d X(> . It can be readily seen that such a conductivity profile 
satisfies the conditions for a steady and quiescent base state. If a* is expressed as 


it follows that 


d 2 o *(z*) 

dz * 2 


= 0 , 


( 14 ) 


and according to Eq. (12) the conductivity profile is steady. 
The base state solution associated with Eq. (13) is: 


E (z ) = 


1 + 


ActY 


"Eq ij ; 




(J 0 )d 


(15) 
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(z 
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1 + 


_J 

A az_ 
ct 0 d 



* * 

Pf(z ) = 



Aov 


} o 


A *A 2 

Act z_ 


(T 0 d 



(16) 


(17) 


Here p‘„ is the ambient pressure. £ 0 is the electric field strength at the lower boundary in 
the base state (from Eq. (15)) and £ 0 a 0 = J 0 , i.e., there is a constant current density J 0 
across the fluid layer in the base state. 

The base-state solution suggests the following reference scales, which are used to 
place the balance laws in dimensionless form: free charge density, £E 0 Aa/da 0 ; electric 
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field, £ 0 ; and length, d. Hereafter, cr denotes the dimensionless conductivity, i.e. 

cr = (a* - <t 0 )/a<t. 

The dimensionless momentum balance reads 

D\ ] - 

~Di = p+ 4g^ v y+PtE - < 19 > 

where the variables are now in scaled form. A balance of the inertial and electrical effects 
requires that the velocity scale be ^eEqAo/pcTq , and the appropriate stress scale is 

eE 0 Aa/a 0 ; the time scale is ^d 2 pa 0 / eEqAcj . The dimensionless group that appears in 
Eq. (19) is defined by 


Gr„ 


d 2 eEq Act 
v 2 p ct 0 ’ 


( 20 ) 


and can be thought of as an electric Grashof number, since Gr e is a measure of the relative 

strength of electrical body forces and viscous effects. The inertial scaling chosen here is 
appropriate when Gr e is 0(1) or more. 

The governing equation for the evolution of the electrical conductivity becomes 


Da 

~Dt 


1 

^eVGre 



( 21 ) 


if an electric Schmidt number is defined as 


Sc c 


v 


(22) 


Here Sc e is roughly 10 3 - 10 6 , since K eff ranges from 1(T 9 - 10 ~ l2 m 2 /s. 
The conservation of charge relation transforms to 


Dpt _ 1 d 2 a 0 

Dt ^Gr c ve | 

where T = d 2 a 0 /ve represents the viscous relaxation time of the liquid normalized to 
e/<J 0 , the charge relaxation time in a liquid of conductivity a 0 and dielectric constant e. 
For electrolytic buffer solutions a 0 = 10 ~ 2 S/m and £ « 1(T 9 C/(mx V), so d 2 oJve is 
O(10 8 ) when d ~ 10 _3 m and v«10 ~ 6 m 2 /s. 

Gauss' law in dimensionless form is 


, Act , 

1 + — <J\p { +E- Vct 


( 23 ) 


V-E 


Ac 



(24) 


Because the electric field is irrotational, the electric potential can be introduced, viz, 


E = ~ V0 , 


(25) 


and then Eq. (23) becomes 



Similarly, substituting Eqs. (24) and (25) into Eq. (19) yields 
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Equations (26) and (27), along with Eqs. (6) and (21), are the balance laws used to 
describe the dynamic behavior of the system whose base state is depicted in Fig. 1. In the 
next section, a linear stability analysis is performed on this base state. Small perturbations 

are introduced to the system and their evolution is analyzed to determine if the disturbances 
grow or decay with time. 

Linear Stability Analysis 

The base-state solution in dimensionless form is: 


<7(z) = z ; 


E(z) = 


1 

i , Act 
1 + — z 


p(z) = 


* 0 / 


, Act - . 
7 : 7T + A) > 


2f , + Acx ^ 
1 + — z 

V ^0 J 


P f(z) = - 


V °o ) 


(28) 

(29) 


(30) 


(31) 


v(z) = 0. 


(32) 
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Note that overbars indicate the base-state variables. 

To cany out a Itnear stability analysis of the base state, the following perturbation 
variables are introduced 


v(x,0 - v'(x,0 - u'(x,t ) + v'(x,t) i + w'(x,t ) i ; 


cr(x,0 = <7{z) + cj'{x,t) ; 


0(x,o = ^(z) + 0 / (x,o ; 


p{x,t) = p(z) + p'(x,t) ■ 


where the prime denotes the perturbation variables, and 


0M = 0o-J o ! E({)rf{ = 0„-^it/ 


, Ac 
1+ — z 

V G 0 y 


Formal substitution of Eqs. (33)-(37) into the balance laws yields 


f=- v+ vsr vV+ S{ gw, l+ gv,l. 


,da 1 , 

— -fw — = — y 2 a - 


* Sc z ^Gr 


— (vV) - d3<l> ■■■' T 


dt 


+ — w — — „ 
dz ^Gr t 


, A <7_ 
1 + cr 

^0 


V 2 0' + 


(33) 

(34) 

(35) 

(36) 

(37) 


(38) 


(39) 
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Act d 2 (j) , + do dty' + d(j> do' 
o 0 dz 2 dz dz dz dz 


V • v' = 0 , 


upon linearization. The pressure and the x and y components of the velocity are eliminated 
from the problem by taking the curl of the momentum balance and making use of 
continuity. The resulting expression is 


~ — V 2 IvV = V 2 (V 2 <p') - V? 0 ' , 

dt ^Gr e J Act dz ' ' Act dz 3 


V 2 m— +il. 

1 dx 2 dy 2 

Equations (39), (40) and (42) are a coupled set of partial differential equations and the 
unknowns to be determined are: v/(x,r), the perturbation in the z-component of the 
velocity; <p'(x,t), the perturbation in the electric potential; and o'(x,t), the perturbation in 
the conductivity. 

Equations (39), (40) and (42) are linear and symmetric in x and y, and, so, are 
amenable to a normal modes analysis. Accordingly the solutions for the perturbation 
variables are assumed to have the form 


vv'(x,t) w'(z) 

(!)'{x,t) > = < <j>\z) >f{x,y)e sl 
o'(x,t)\ <7'(z) 
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There are solutions of similar form for p'(x,t), u '(x,t ), v'(x,t) and E'(x,t). Formal 
substitution of Eq. (43) into Eqs. (39), (40) and (42) yields 



Op d<P 
Act dz 




A<t dz 3 



(44) 





(45) 



T 

V^e 


, Acr_ 
1 + <7 

°0 


t Act 



da d[b'fe st ) 
dz dz 


+ 




dz dz 


(46) 


Based on Eq. (45), a form of the disturbances in the x and y directions can be 
determined 19 , viz. 


/(*.*) *«' (B ' X+<V) 


with 


a 1 = a\ + a* . 


Equations (44)-(46) then become, 


* - - ° 2 ) r - - -« 2 S # (° 2 - « 2 )^- + • 


Ac cfe 


Act dr 


(47) 
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, dz 
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From the base state solutions, 


d<J 

dz 


= 1 ; 
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i Act 
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Substitution of the above relations into Eqs. (47)-(49) yields 
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( 51 ) 
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The above set of equations, Eqs. (50)-(52), is a coupled system of ordinary 
differential equations in z; two equations of the set have coefficients that vary with z . The 
set of ODEs is solved with the following choice of boundaiy conditions: the electric field, 
velocity and conductivity are fixed at their base-state values on the upper and lower 
boundaries, i.e. dty'/dz = w' = o' = 0 at z = 0 and 1 , and the boundaries are either stress- 

free, i.e. d 2 w'/dz 2 = 0, or rigid, i.e. dw'/dz = 0. These boundary conditions are 
consistent with the constant current imposed in the base state: J* (z *) = J 0 = E 0 a 0 for 

0 < z <d. Note that J' = 0 at z = 0 and 1 . The eigensolutions to Eqs. (50)-(52) may be 
determined numerically for each wave number a, but some insight into the stability 
problem may be gained by manipulating the equations further. 

First let r -> oo, which is a reasonable approximation because the relaxation time of 
the charge e/a Q is of 0(1 0 8 ) times faster than the characteristic relaxation time of the 
fluid. Then, for neutral stability (s s 0), Eqs. (50)-(52) read 




1 A(J 
1 + — z 

^0 
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Substituting Eq. (54) into (53) yields 
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(° 2 - 2 ) 3 ^ 


ArJ f Ac 7 ' 
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(z) 2 -a 2 )0'-2a : 


Act 


1 


°o [ Act | 

r*‘j 


73 0'- (56) 


Notice Eqs. (55) and (56) can be solved for o' and p' without reference to vv'; the only 

parameters that appear in these two equations are Act/ct 0 , Sc e Gr e and cc. Moreover Eq. 

(54) shows that in order for w' * 0, it must be that o'* 0 . Thus the conditions for 
marginal stability (with a finite and t-+oo) depend only on three parameters, viz. 
Ao/o 0 , Sc e Gr e and a. This is confirmed by the numerical results shown in section VI. 

The product Sc e Gr e can be thought of as an electric Rayleigh number, 
Ra e = Sc e Gr e = d 2 eEqAo } [iK^Oq . If the characteristic velocity ^eEqAo/oqP is 
denoted by u c , then one interpretation of Ra c is as follows: (eE^Ao/o 0 )/(pu c /d) is the 

ratio of electrical to viscous stresses; and u c d/K eff is a Peclet number that represents the 
relative strength of convective and diffusive transport of the conductivity fluctuations. So 

Ra = g£ 0 Acr/gf) u c d 
C f*c/d K t ff 



is a parameter that weighs the destabilizing effects of the electrical stresses against the 
stabilizing influence of viscosity, as well as the mechanisms that govern the evolution and 
transport of conductivity perturbations: convection (destabilizing) and diffusion 
(stabilizing). 

Numerical Approach to the Neutral Stability Analysis 

To generate the neutral stability curves for Eqs. (50)-(52), the conditions for which there 
are eigensolutions to Eqs. (50)-(52) (with j = 0) were found by the following numerical 
technique. Equations (50)-(52) were first re-cast as a set of eight first order ODEs in z. 
The general numerical solution to the system of eight ODEs was constructed for specified 
values of Sc e , Gr c , Act/cTq and a. This was done by solving the system of equations for 

each of eight linearly independent sets of initial conditions (prescribed at the lower 
boundary, z = 0) and the integration was carried out using the IMSL Library routine 
DIVPAG. If y- denotes the / th linearly independent solution to the system of eight ODEs, 

the general (numerical) solution to the system can be written as 

8 

y=5>/y«’ 

i=l 


where the c i are arbitrary constants. The general solution y must satisfy, for example, the 
following set of conditions at the (stress-free) boundaries: 
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dp 

dz 


'z=0 


_ dp 

dz 


= 0 


z= 1 


w{z = 0) = w'(z = 1) = 0 
&'(z = 0) = &'(z = 1) = 0 


d 2 w' 

~dF 


z=0 


d 2 w' 

dz 2 


= 0 


z=l 


(57) 


The task of satisfying Eq. (57) consists, at least conceptually, of choosing the appropriate 
values for the c f (/ = 1,2,...,8) . To determine these values, one must solve an 8x8 

system of linear algebraic equations, which can be expressed in the form 


Ac = 0, c = 


The elements of the matrix A depend on y, (0) and y,(l), the linearly independent 
solutions evaluated at z = 0 and 1, respectively, as well as Sc e , Gr c , etc. Now the ODEs 

and the initial conditions are a set of homogeneous equations. As such, solutions for 
& (z), ^ (z) and w (z) that satisfy Eqs. (50)-(52) and Eq. (57) are generally trivial, i.e. 


Ac = 0 => c = 0 => y = 0. 

However for certain combinations of Sc e , Gr e , etc., nontrivial solutions do exist; these are 
commonly known as the eigensolutions to the problem and they obtain when det A = 0 . A 



search of the parameter space was implemented to establish the conditions for which 
det A = 0 ; this was done by varying Gr e , while holding all other parameters fixed. The 

conditions thus delineated are used to plot the neutral or marginal stability curves presented 
in the following section. 

Results and Discussion 

Results for the neutral stability curves are shown in the figures and tables that follow. The 
marginal stability curves for a fluid layer with Aa/cT 0 = 10 and stress-free boundaries are 
shown in Figs. 2 and 3. In Fig. 2, Gr e is plotted versus a at Sc e =10 3 , 10 4 and 10 5 . 
In Fig. 3, these data are re-plotted as Ra c versus a, and the three curves collapse into one. 

Thus, as indicated in section IV, the marginal stability curves associated with Eqs. (50)- 
(52) depend on Acr/cr 0 and Sc t Gr t when T -> «>. These results are also shown in Table 
I, where marginal stability values of Ra e vs. a are listed for Sc e = 10 3 , 10 4 and 10 5 ; 
Ra e values agree to four digits across a two decade variation in Sc c . For A<j/<j 0 = 10 and 
stress-free boundaries, the critical Ra t is 1.504 x 10 4 and the critical wave number is 1.98. 
When A<r/cr 0 is increased, the minimum of the marginal stability curve shifts to 

slightly higher wave numbers (Fig. 4). More importantly, the magnitude of the critical 
Ra e increases with Acr/<J 0 . This seems counterintuitive inasmuch as the region of 

stability widens as the conductivity increment of the base state increases. Recall, however, 
that 

°0 

So, even though the critical Ra c increases monotonically with A<x/ct 0 , the critical electric 
field strength does not. Plots of T 2 = eE 2 d 2 /nK t({ versus a are shown in Fig. 5; T 2 can 
be thought of as the square of the dimensionless field strength or, alternatively, as the 


scaled electric energy density. Notice that the magnitude of the critical £ 2 goes through a 
minimum as Act/cTq is increased. This again seems to contradict the notion that 

sharpening the conductivity increment in the base state favors instability. However, it will 
be seen that the widening of the stability region for Ao/cr 0 > 20 is due to boundary effects. 

The imposition of the electric field on the base state conductivity gradient induces 
within the fluid layer a volume charge, the distribution of which is given by Eq. (31). The 
base state field acts on the volume charge to produce a coulombic body force that is 
proportional to E z dE z /dz. This electrical body force is opposed by the base state pressure 
gradient. But, for reasons outlined in Section I, this balance of forces is not necessarily 
stable. When a small fluid parcel is displaced upward (or downward) it is surrounded 
locally by a fluid of higher (or lower) conductivity. Due to the conductivity differences, a 
polarization charge, proportional to the local strength of the base state electric field E , is 
induced about the periphery of the parcel. This polarization charge develops 
instantaneously compared to the response times of the fluid and the conductivity profile, 
i.e. 

e/<J 0 « d 2 /v and d 2 /K eff . 


The local field gradient acts on the polarization charge much like it would on an induced 
dipole, exerting a force that is proportional to E z dE z /dz. This force is similar to that 
which gives rise to the dielectrophoresis of leaky dielectric drops and tends to perpetuate 
the motion of the parcel in the direction of its original displacement^. If the viscous 
stresses can resist the motion of the parcel well enough for ionic diffusion to dissipate the 
local conductivity differences, the motion of the parcel is arrested. If they cannot, the parcel 
continues on upward (or downward) and instability ensues. 



In Fig. 6, the profile of E z dEjdz is plotted for several values of A cr/<r 0 . The 
portion of fluid in which the (destabilizing) dielectrophoretic forces are appreciable is 
proximal to the lower-conductivity boundary (z = 0) and becomes increasingly narrow as 
A<J/<J 0 increases. Accordingly, at high values of A<x/ct 0 , one would expect the lower 

boundary to exert a strong stabilizing influence, especially in comparison to that of the 
upper boundary. 

Figures 7, 8 and 9 are plots of the marginal stability curves £ 2 vs. a at 
Ao/( 7 0 = 1000, 100 and 50, respectively. Results are shown for four kinds of boundary 

conditions: (1), both boundaries rigid; (2), the upper boundary stress-free and the lower 
one rigid; (3), the upper boundary rigid and the lower one stress-free; and, (4), both 
boundaries stress-free. The effect of the conditions imposed at the upper boundary is 
clearly secondary to that of the conditions imposed at the lower. This is evident in two 
ways. First, the marginal stability curves of those cases that share a common boundary 
condition at z = 0 are closest. Second, changing the boundary condition at z = 1 (e.g. 
from stress-free to rigid) shifts the marginal stability curves by a factor of 2 or so in £ 2 , 
while a change in the condition at z = 0 shifts the curves by as much as a factor of 100. It 
would seem then that the minimum in the variation of £ 2 rit with Aa/o 0 (Fig. 10) is due to 
the increasing influence of the lower-conductivity boundary as Aa/a 0 increases. 

Comparison of the foregoing results with experimental observations are favorable 
but far from definitive. Electrohydrodynamic instabilities were manifested in microgravity 
IEF experiments 10 performed on aqueous electrolyte buffers at field strengths of 
approximately 1 .IkV/m and current densities of 17A/m 2 . Numerical simulations of the 
IEF processes indicate that, were these IEF experiments to have proceed without the onset 
of convection, conductivity variations exceeding an order of magnitude would occur over 
distances of about 1 mm 11,12 ’ 15 . Using the viscosity and dielectric constant of water, 



d = \mm, ^ = 10 10 m 2 /s, and a conductivity <r 0 = 1(T 2 S/m, the results for stress-free 
boundaries and Ac 7 / cr 0 = 10 give a critical field strength £ 0 crit = OAkV/m and a critical 
current density J 0 crit =4A/m 2 . Thus, under experimental conditions for which the 

predicted critical field strength (or, alternatively the critical current density) was exceeded, 
instability was indeed observed 10 . 

Concluding Remarks 

This analysis shows that electrically-driven convection can ensue in an initially quiescent 
fluid layer with an electrical conductivity gradient. The fluid instability is due to 
dielectrophoretic forces that stem from local electrical conductivity differences in the fluid. 
If diffusive transport of charge-carrying solutes is sufficiently fast, the dielectrophoretic 
forces vanish with the relaxation of conductivity differences and quiescence persists in the 
layer. In this work, the diffusive relaxation of conductivity differences has been lumped 
into a single parameter K e{{ , which appears in a convective-diffusion balance on < 7 . 

The implications of the results are that, in order to maintain stability in the fluid 
layer, the magnitude of the conductivity gradient and electric field must be controlled. The 
magnitude of A<t/c 7 0 dictates the magnitude of the critical field strength that can be 

imposed across the layer before electrohydrodynamic instability ensues. Moreover, it is 
shown that the conditions prescribed at the lower-conductivity boundary have a significant 
effect on the stability of the fluid layer. 
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Table I 

Marginal Stability Values for Ra e versus a and Sc e at Act/<t 0 = 10. 


a 

Ra c @Sc e = 10 3 

Ra e @ Sc e = 1 0 4 

Ra c @ Sc e = 1 0 5 

1.0 

25460 

25460 

25463 

1.1 

22512 

22511 

22508 

1.2 

20339 

20339 

20344 

1.3 

18721 

18721 

18718 

1.4 

17513 

17513 

17510 

1.5 

16618 

16617 

16619 

1.6 

15967 

15968 

15971 

1.7 

15515 

15514 

15510 

1.8 

15226 

15226 

15224 

1.9 

15077 

15076 

15081 

2.0 

15048 

15049 

15048 

2.1 

15128 

15129 

15125 

2.2 

15308 

15307 

15312 

2.3 

15579 

15579 

15576 

2.4 

15939 

15939 

15938 

2.5 

16385 

16385 

16389 




Fig. 1 An initially quiescent and steady fluid layer of infinite extent in the x and y 
directions with a linear electrical conductivity gradient in the z direction. 
Asterisks are affixed to the variables to indicate that they are in dimensional 
form; cf. Eq. (13). 



Transverse wave number, a 


Neutral stability curves for stress-free boundary conditions; Acr/c t 0 = 10; 
t- 10 . Legend: 1. Sc e = 1 0 3 , Gr c crjt = 15.04, a crit = 1.98; 2. Sc e = 10 4 , 
Gr c. crit = 1-504, « crit = 198; 3. Sc c = 10 5 , Gr e cri( = 0.1504, a • =1.98. 





Transverse wave number, a 


Neutral stability curves for the stress-free boundary conditions; r = 10 8 . 
Legend: 1. Ao/<t 0 = 10, Ra t crit = 1.504 x 10 4 , a crit = 1.98;2. Acx/cr 0 = 20, 
^e.crit = 2.430 xlO 4 , a crit =2.06; 3. Acr/<r 0 =50, Ra e crit = 6.788 x 10 4 , 
^cnt =2.16; 4. A<7/(7 0 = 100, Ra t crit = 1.766 x 10 5 , « crit =2.22; 5. 

Act/ct 0 = 1000, Ra. crit = 7.407 xlO 6 , a rril =2.34. 



Transverse wave number, oc 


Neutral stability curves for the stress-free boundary conditions; r=10 8 . 
Legend: 1. A<r/<r 0 = 10, £ 2 rit = 1.505 x 10 3 , a crit =1.98; 2. Act/(t 0 =20, 
■E crit = 1.215X10 3 , a crit = 2.06; 3. Act/ct 0 = 50, £ 2 nt = 1.358 x 10 3 , 
«cm=2.16;4. A<t/(7 0 = 100, £ 2 rit = 1.766 x 10 3 , a crit =2.22;5. 

A<?/<t 0 = 1000, £ 2 nt =7.407xl0 3 , •, =2.34. 
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Fig. 7 Neutral stability curves for A<t/<t 0 =1000 for the four kinds of boundary 
conditions. Legend: 1. rigid-boundaries; 2. the upper boundary stress-free and 
the lower one rigid; 3. the upper boundary rigid and the lower one stress-free; 
4. stress-free boundaries. For all these curves t = 10^ 
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Transverse wave number, a 


Fig. 8 Neutral stability curves for Act/< 7 0 = 100 for the four kinds of boundary 
conditions. Legend: 1. rigid-boundaries; 2. the upper boundary stress-free and 
the lower one rigid; 3. the upper boundary rigid and the lower one stress-free; 
4. stress-free boundaries. For all these curves r = 10 8 . 
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Transverse wave number, a 


Neutral stability curves for Ao/<r 0 =50 for the four kinds of boundary 
conditions. Legend: 1. rigid-boundaries; 2. the upper boundary stress-free and 
the lower one rigid, 3. the upper boundary rigid and the lower one stress-free; 
4. stress-free boundaries. For all these curves t= 10 8 . 





Appendix 1 : An Estimate of the Conductivity Gradients 

Required to Produce the Convection Observed 
in IEF Space-Flight Experiments 

To obtain an estimate of the magnitude of the conductivity gradients needed 
to yield the flow strengths observed during the (STS- 11 and STS-26) 
microgravity experiments, a scale analysis of the momentum balance is 
performed as follows. The Reynolds number is a measure of the relative 
magnitude of inertial and viscous effects, and it is defined as 

■ ■ , p ^v 

Re = inertial = /T _ _ U c l c 

viscous jj,U c / l*-/ ~ v (i) 

M /p 


For the present analysis, the focusing chamber is assumed to have 
cylindrical geometry, and from experimental results, the flows encountered 
have a velocity on the order of 1 cm/min, in a capillary of 0.5 cm 
diameter. The kinematic viscosity of the buffer solution is approximately 
that of water at room temperature, lO^mV. Thus, 


Re _ (10~ 2 /«/min)(5x1Q~ 3 m)(^min/.y) 

10 ~ 6 m 2 /s 


(2) 


which indicates that the inertial and viscous effects are of comparable 
magnitude. 


To asses the significance of the electrical effects, they are compared 
to the viscous effects, since the momentum generated by electrical stresses 
is transferred at the onset of motion by viscous action. From the Navier- 
Stokes equations, a balance of the viscous and electrical terms requires 

/xV 2 v -Vp~ ^-E • EVe-p f E ( 3) 


where « is understood to mean that the terms are of similar magnitude. 
Substituting Eqs. (8) and (11) from chapter 2 into (3) above, one finds 


< uV 2 v - V/7 = ^-E • E Ve • 


^v-V(Ve-E + eV-E) + ^Ve-e~J.E 


E. (4) 


Since for IEF the gradient in the dielectric constant is negligible, 


fiV 2 v-Vp 


£ 

— V 
(7 


•V(eV- E)-e~ 

<j 


E 


E. 


(5) 


Also 


V(eV • E) = Ve(V • E) + eV(V • E) s eV(V • E), 


( 6 ) 


which leaves 


From Eq. (7), it follows that 


Vo 


U 


£ ~K~n^ + -u^- 


( 8 ) 


or 


L + i!L, 

ct M Z 2 e£ 2 o Z 2 


(9) 


For the IEF experiments ^a/(eE c f » l, so 


Vo 


(10 3 fcg/(mx,y))( 10 2 m/min)(^min/^) 

C /I w 1 A“I2 //l/i, . \\/i /'/•n ¥ 7 I \ 2 . - 


cr (80)(8.854 x 10“ 12 C/( V x m))(l667 V/m) 2 (5 x 10" 3 m) 2 


— 3 m 


-i 


( 10 ) 


Equation (10) indicates that Vo/o need only be ca. 3 m' 1 to produce flows 

of strength 1 cm/min in the IEF flight experiments. The simulation results 

shown n Fig. 21 of Chapter 1 indicate that gradients of this magnitude are 
realizable in IEF. 
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Appendix 2 : The Charge Conservation and Conductivity 

Evolution Equations 

Conservation of Charge 

In systems where conduction involves more than one charged species (e.g. 

aqueous electrolytes), each of the charges is governed by a conservation 
equation 


-^ + V-(p ± v ± J ± ) = 0, (1) 

where the current density for a bipolar conductor is 

J±=b ± p ± E + K ± Vp ± , (2) 

Here b ± are the mobilities and K+ the diffusivities of the cation (+) and 
anion (-). The net free charge density and conductivity are 


Pt-P+-P- ; ° = KP+ + h -P- , 


(3) 


and 


P± = 


b + +b_ 


(4) 


Substitution of Eqs. (2) and (3) into Eq. (1) and introducing Gauss' law, 


viz. 



